Single Polymer Confinement in a Slit: Correlation between structure and dynamics 
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In this paper, we construct an effective model for the dynamics of an excluded-volume chain under 
confinement, by extending the formalism of Rouse modes. We make specific predictions about the 
behavior of the modes for a single polymer confined to a slit. The results are tested against Monte 
Carlo simulations using the bond-fluctuation algorithm which uses a lattice representation of the 
polymer chain with excluded-volume effects. 



INTRODUCTION 

Confining polymers inevitably leads to a competition 
between length scales arising from internal cooperativ- 
ity and those imposed by external geometry, resulting in 
qualitative changes of the static and dynamic behavior. 
The effects of confinement on a single polymer chain is 
still not completely understood. A scaling picture has 
been used to describe the changes in the relaxational dy- 
namics of a single, self-avoiding chain confined to a tube 
or a poreP, Since a self- avoiding chain under cylin- 
drical confinement models a variety of phenomena rang- 
ing from DNA translocation through synthetic pores [3| 
to bacterial chromosome segregation 3, there have been 
a number of simulation studies investigating its statics 
and dynamics [E S S 0, S i] • 

In this paper, we use the- 
ory and simulations to study the dynamics of a single, 
self- avoiding polymer confined to a slit. 

One of the key concepts in polymer dynamics is the 
relaxation of modes that describe the internal dynamics 
of a chain molecule [3, [13, 11 1. We construct a theoretical 
framework for describing the relaxation of the modes of 
a confined, self-avoiding chain in the absence of hydro- 
dynamic interactions, by extending the Rouse model [lOt- 
Our analysis shows that the anisotropy, introduced by 
the slit geometry, leads to dramatically different relax- 
ations of the longitudinal and transverse modes. This 
effect, in turn, changes the nature of the motion of sin- 
gle monomers, and the collective dynamics of the self- 
avoiding chain. We perform simulations using the bond- 
fluctuation model[l^, and show that the results can be 
explained, semi-quantitatively, by our theoretical model. 
Explicit theoretical results for the relaxation of Rouse 
modes, presented this work, provide a powerful tool for 
predicting fine scale and large scale motion of biopoly- 
mers within the confines of a cell. 



ROUSE MODES 

If excluded volume and hydrodynamic interactions are 
ignored, the dynamics of an unconfined, ideal chain 
(without self-avoidance) is well described by the Rouse 



model0,[ll: 



(-Qf ^ fc(Xn+l - 2x„ 



x„-i)+g„(t) (1) 



where x„ denotes the position of the n monomer of 
size h, k — , g„ {t) is a random force arising from the 
solvent surrounding the polymer, and C is the friction 
constant of a single monomer. Considering dynamics at 
length scales much larger than the individual monomer 
size b. 



and Eq. ([T]) becomes a partial differential equation: 



c 



dx{n,t) d'^x{n,t) 



dt 



(2) 



(3) 



The random force, g{n,t), is Gaussian and is character- 
ized by the following correlations: 



< gjin,t) >= 



(4) 



< 9j{n',t)gk{n',t') >= 



N 



S,kS{t~t')S{n-n') (5) 



with j and k being the cartesian components of g(n,t), 
and N the number of monomers in the chain. The normal 
modes of the above equation are the Rouse modes, and 
they decay exponentially with characteristic time scales 
that depend on the mode index 13, 13] ■ 

A major missing piece in the above description is ex- 
cluded volume interactions [l3| • The inclusion of these in- 
teractions (or confinement) introduces couplings between 
the Rouse modes, and is described by an equation of the 
form: 



c 



d-x.{n,t) d'^x{n,t) 



dt 



-g(x,i)+F(n,t) (6) 



where F(n,<) represents the force on the n monomer 
due to all the other monomers, or due to the confining 
geometry. The presence of this term couples the Rouse 
modes, and in general, this equation has no analytic so- 
lution. An approach that has been used for unconfined. 
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self- avoiding, polymers is to assume that the dynamics 
can still be described by decoupled Rouse modes [l3|, but 
with renormalized parameters. Bond fluctuation algo- 
rithm simulations indicate that for such polymers, the 
Rouse modes are indeed weakly coupled [1^] and that 
the renormalized Rouse modes describe the dynamical 
behavior. The relaxation of the modes is characterized 
by a time scale Tp ~ ^„+i , where the mode index p = 

0, 1, 2, , and v is the Flory exponent: u = 3/{d + 2) 

with d being the spatial dimension [lo|. In 2D, Tp ~ ^^kjt- 
The longest relaxation time is ri oc N'^'^'^^, and is associ- 
ated with the reorientation or rotation of the polymer [loj. 

RelcLxation of Rouse Modes for Confined Polymers 



The generalized Rouse model treats a confined, self- 
avoiding polymer as an ideal polymer with renormalized 
spring constants. The utility of the approach relies on ac- 
cess to reliable results for the static correlations between 
monomers, < {xmj — Xnj)'^ >eq- The time scales associ- 
ated with the relaxation of the Rouse modes, and ulti- 
mately the dynamics of the chain and single monomers 
can be obtained, once these static correlations are known. 
The task at hand is, therefore, to find an appropriate 
model for the equilibrium correlations for confined poly- 
mers and to predict the dynamical behavior using Eqs. 

©-(HBO. 

The next section is devoted to discussing appropriate 
forms for static correlations function of a confined, self- 
avoiding polymer. 



In this work, we construct a framework for describing 
Rouse modes in confined polymers. To facilitate discus- 
sion, we briefly outline the steps leading to the results for 
unconfined polymers with self-avoidance. Neglecting the 
coupling between Rouse modes leads to a set of linear 



equations for the mode amplitudes Xpj 13, l3 |: 



' dt 



kpjXpj + gpj{t) 



(7) 



Here 



Xr,i — 



N 



N 



dn < 



'\^)^n,{t), (8) 



and (p — 2N( for p > 0, — N^. All the self-avoidance 
and boundary effects are incorporated through the kpj 
term. Requiring that in the long time limit, the Rouse 
modes have the correct equilibrium distribution, implies 
thatfH: 



< ^pj >eq 



(9) 



where < . >eq is an equilibrium average. The exact ex- 
pression for < Xp.j >eq can be derived using Eq.®: 



(^) ^ < {Xno-X^j f > (10) 



As shown by Doi and Edwards |13l|. for p » 1, < Xpj >eq 
can be expressed as: 



< > 



N 



P3 ^^-3 ^p2^2 



N 



ducos y-j;^ J X 



d <^ i^Xraj •^Tij') ^ eq 

dndm 



(11) 



where u = m — n. 



BLOB MODEL FOR SLIT GEOMETRY 

For an unconfined polymer, consisting of N monomers 
of size 6, represented as a self-avoiding- walk (SAW), the 
radius of gyration scales asflij: 



Rg - bN" 



d+2 



(12) 



Upon confining this SAW to a slit, we expect the ex- 




FIG. 1: Sample configuration from simulation of a 2-D poly- 
mer confined to a tube. Circles indicate rough shape of the 
blobs and have a size on the order of the tube width, L. 

eluded volume effects to change the monomer-monomer 
correlation functions. The scaling of these correlations 
with the tube width and length of the polymer has of- 
ten been described by the blob modeljj]. Mean-field 
approaches have also been used to model these static 
correlations [1, [l3| . 

In a slit of width L, the first effect of confinement is 
to make the polymer conformations anisotropic [l^, and 
the shape needs to be characterized by a transverse {R±) 
and a longitudinal (i?||) length scale. For widths, L < Rg, 
where Rg is the radius of gyration of the unconfined poly- 
mer, confinement has the effect of stretching the polymer 
in the axial direction and screening out the excluded vol- 
ume interactions beyond a length scale comparable to 
L. The polymer can be visualized as a succession of 
blobs with excluded-volume effects being maintained for 
distances smaller than the blob sizeQ. The number of 
monomers per blob, g, can be calculated by requiring 
that the radius of gyration of these monomers matches 
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the blob length scale [3 For the slit geometry, the 
length scale of the blob is assumed to be L|15I| , and there- 
fore, 



The polymer is composed of iVt, 
lateral extension of: 



(13) 

N / g blobs and has a 



R\\=LNb = Nbib/Lf'"-^ ^ Nb{b/Lf/^ , (14) 

where the last relation is specific for the slit geometry 
(2D). Since the tube width confines the polymer in the 
perpendicular direction[21]: 



R±(xL 



(15) 



In papers by Edwards and Singh, estimates of polymer 
sizes were made using a mean- field approach 17|. This 
method was later applied to a polymer in a tube by Mor- 
rison and ThirumalaiQ. The blob predictions for the size 
of a polymer in narrow tubes agree with their estimates 
up to a factor of order unity. 



Monomer-Monomer correlations in the Blob Model 

The arguments presented above hold equally well for 
sub-chains as long as the sub-chains we consider are 
large enough to have the asymptotic behavior of a SAW. 
Therefore, the correlations between monomers in the di- 
rections parallel and perpendicular to the slit walls are 
should scale as: 



\n 
In 



mp/^fo^ |n - m| < ( j)^^^ 

m|268/37,~2/3 \n-m\> (^)^^' 
(16) 



13/2^2 



4/3 

(17) 

Recent simulations [15| have measured the distribution 
of blob sizes by analyzing sub chains, and for strong con- 
finement {L < Rg), their results are in agreement with 
the blob model predictions. In order to proceed with the 
dynamical calculations, we assume that the correlations 
are piecewise continuous with sharp transitions between 
the different scaling regimes represented in Eqs. (jl6p and 



Simulation results for the Monomer-Monomer 
correlation function 

We have used simulations to measure the monomer 
correlations and compared them to the predictions of the 
blob model outlined above. To our knowledge, these pre- 
dictions have never been directly confirmed. 



t 



FIG. 2; Example chain configuration. Monomers occupy each 
box. Allowed moves are labeled with dashed arrows. Bonds 
are represented by solid lines. 



Simulation Technique The simulations were based on 
the bond fluctuation algorithm, a course-grained lattice- 
based algorithm that allows for the analysis of dynamical 
properties of polymers T^. As the name implies, the al- 
gorithm allows for fluctuations of bond lengths and em- 
ploys only local moves of monomers. Excluded- volume 
effects arc taken into account by ensuring the initial poly- 
mer configuration contains no intersecting bonds and by 
allowing only the following set of bonds: 

S = (2,0)U(2,1)U(2,2)U(3,0)U(3, 1)U(3,2) 

This also leads to an effective monomer size 6, where 
b « 2.8 lattice spacings. This algorithm has been 
shown to reproduce the right scaling for polymer size 
and follow the generalized Rouse dynamics for excluded- 
volume polvmers[lo| in all spatial dimensions. It has 
been widely used to study the dynamics in different 
environments [lil 22, 23, 24 1. In this paper, the simu- 
lation results are used to test the predictions of our the- 
oretical formalism, and to explore in detail the dynamics 
of a chain confined to a slit. 

Results Monomer correlations were obtained for 
chains with N = 200, Rg = 149. This size is optimal 
in the sense that it is large enough to illustrate the ef- 
fects of blobs while still yielding relaxation times that can 
be accommodated within accessible simulation times. As 
will be shown from our calculations of the Rouse mode 
relaxations, the longest time scale increases as and 
as L~^/^, making it difficult to determine correlations 
functions for long chains in a narrow slit. 

The initial chains were generated using a back-tracking 
algorithm which would redraw the part of the walk if a 
new segment had not been placed within 30 attempts. 
To ensure little biasing from the generation of the initial 
walk, chains were equilibrated by allowing the chain to 
undergo 10^ monte carlo steps (mcs) (200*10^ attempted 
moves). After equilibration, data was sampled every 10"^ 
mcs for runs of 10^ mcs. The process was repeated for 10 
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initial configurations and averaging was conducted over 
all sampled configurations of all runs as well as over the 
whole chain. This averaging leads to the data shown in 
Figs. Handgl 

The scaling laws observed in simulations are in broad 
agreement with the blob picture for the monomer cor- 
relations perpendicular to (Fig[3]) and parallel to (Fig|31) 
the walls. In reporting the results of the simulations, we 
measure all lengths in units of b. Fig. [3] demonstrates 
that the data for the transverse correlations for different 
values of L can be made to scale if |n — r7i| is divided by 
L'*/'^ and the correlations are scaled by L^, confirming 
the form of Eq. \T7\ and the scaling of the blob size g{L) 
as L^/'' (i/ = 3/4 in 2D). The results for the correlations 
in the longitudinal direction indicate a crossover between 
two power laws, as shown in Fig. |4l 

According to the authors' knowledge, this is the first 
direct test through simulations of distinct scaling regimes 
predicted by the blob model. Since there is broad sup- 
port for the blob-scaling picture of monomer correlations 
from the simulations, we proceed to utilize these results 
for calculating the properties of the Rouse modes under 
confinement. 
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FIG. 4: < {x„\\ — K^ii)^ > versus slit width for the same set 
of parameters as Fig. |3] A crossover between two power laws 
is observed. 




|n - m|/L* 



FIG. 5: Scaling of the blob size, g{L), extracted from the 
crossover points of < {x„\\ — Xm\\)^ > shown in Fig|4] The 
best fit is for L^, however, the L'^^^ behavior evident in the 
transverse correlations is within the error bars on this data. 



FIG. 3: Scaling of < {xr, 



_) > , demonstrating the the 



scaling of g(L) as L"^^^, and validity of Eq. 1171 The x-axis 
has been scaled by L'*''^ , and th y-axis has been scaled by 



Eq. (jTO|) . For these modes, 
where the constant A = 



> 



The chain length is N 
are 20, 30, 40, and 50. 
the lattice spacing 



= 200, Rg 

All lengths are measured in units of 



eq— A-B{p/pmax) , 

In the regime where 
satisfy the condition 
for the use of Eq. pT|) . since these mode numbers satisfy 
p > N{h/Rg f/^ = N. Combining Eqs. ^ and ^ 
we, therefore obtain, for p > Pmax, 



149, and the slit widths 1^ « Rg, modes with p > p„ 



RELAXATIONS OF ROUSE MODES 



< Xpj_ >eq- 



N 



Ap'^TT^ 



(L/b) 



4/3 



du cos 



/pTTU\ 



V N J uV2 



(18) 



To study the dynamical correlations in the anisotropic 
geometry of the slit, we analyze the transverse modes, 
Xp± , and longitudinal modes , Xp^ | , separately. 

Transverse Modes The transverse modes correspond 
to taking j to be in the direction perpendicular to the slit 
axis. The blob picture provides a natural scale of Pmax = 
N/g = Nb'^^^/L^^^, separating the large p from the small 
p regime. For p « Pmax, we use Eq. (fT7|) directly in 



In the above equation, and below, we set 6=1, and 
C/ksT — 1. The relaxation of the transverse modes is 
exponential with a time scale: 

Tp± = ^^N< Xp\ >eg (20) 
Kp± 
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A reasonable approximation of Fc{x), the Fresnel Cosine 
function, is a linear function of the argument for x <^ 
1 and 1/2 for x ^ 1. For p > Pm(ja;,the argument of 
Fc{x) is always greater than 1, and, therefore, Tp± scales 
as which is a form identical to the relaxation of 

unconfined modes. This result provides a self-consistency 
check on our calculations, since the blob model assumes 
that the correlations within the blob are identical to that 
of an unconfined chain. 

These results show that < X^^ >eq decreases mono- 
tonically with p, initially as A — Bp^, and then as p^/^. 
The slowest transverse mode is, therefore, p — 1, and the 
associated relaxation time is given by : 



(21) 



Longitudinal Modes For modes with p > Pmaxi the 
functional form of the mode amplitudes is the same as 
that for the transverse modes: 



<^Pll> 



23/2^5/27^2 



2p 



P-n 



(22) 



For p < pmax , we use Eq. ^61) directly in Eq. (|T0| . In the 
limit of small slit widths, L ^ Nb, we obtain a closed 
form expression [25}: 



'eq' 



(1 - i-l)P)N^ 



(0)L 



2/3 



(23) 



In the small p regime, therefore, the odd and even modes 
behave differently. The odd modes exhibit a p^ decay of 
< X^u >eq, and dominate for small L, since the leading 



VII 
term in Eq 

Eq, 



^vanishes for the even modes. For p > N/ g, 



imphes that < 



> 



eq 



decays as p^/^ 



The longest relaxation time in the longitudinal direc- 
tion is determine by the relaxation of the p = I mode, 
and is given by: 



rll = 2N^/L 



2/3 



(24) 



Unlike t^, this time-scale increases as the slit width L 
decreases. The equilibration time is dominated by the 
relaxation of the p ^ 1, longitudinal mode, and as we will 
show from the simulations, corresponds to a collective 
rearrangement of the chain where the end-to-end vector 
changes sign. 

In the next sections, we will use the results for the 
Rouse modes to analyze the dynamical properties of the 
confined chain. 



Collective Dynamics 

The longest relaxation time of an unconfined polymer 
chain is associated with the relaxation of the end-to-end 
vector, P(t) — RN(t) — Ro(t) and is reflected in the 
autocorrelation function: 



C(t) = (P(t).P(0)) 



(25) 



The Rouse model result for an unconfined, self avoiding 
polymer is[loj: 



C{t) ~cxp ( - — 



(26) 



where ti is the relaxation time of the p ~ \, unconfined 
Rouse mode: ti ^ A^3/2 2D[i3|. In a slit, we need 




FIG. 6: Simulation Results: End-to-end correlation function 
in the perpendicular direction for a chain 50 monomers long 
and for slit widths of 20, 30, 40, and 5000 lattice spaces. 

to look at the relaxation of the perpendicular and par- 
allel components of P(t), separately. Writing C_i(<) and 
C||(i)) in terms of the Rouse modes [lo|, it becomes evi- 
dent that the long-time relaxation is exponential with a 
time scale given by the longest-lived transverse and longi- 
tudinal modes, respectively. From the analysis of Rouse 
modes for the confined polymers, it follows, therefore, 
that 



CAt) 
C\\{t) 



exp 1 --^ 



exp 



t 



= exp 



exp 



110/3 
2m 



(27) 



These results demonstrate that the collective dynamics of 
a polymer confined to a slit is qualitatively different from 
an unconfined polymer. The relaxation time of the lon- 
gitudinal component of the end-to-end vector increases 
dramatically as L is decreased, diverging as L ^ 0. This 
increase in relaxation time has been observed in earlier 
simulations [isj and, as shown in Figs. [7]and[5]is captured 
by our simulations. Since the time scale associated with 
changes in orientation of the end-to-end vector increases 
as N'^L^'^/^, long Monte Carlo runs are required to en- 
sure that chain configurations and their mirror images 
are sampled equally. In order to keep the simulations 
manageable, we studied relatively short chains, = 50 
monomers long. The results were averaged over 10^ mcs 
after initial equilibration of 10^ mcs steps and for 5 dif- 
ferent initial configurations. 

In the unconfined polymer, the slowest mode corre- 
sponds to the reorientation of the end-to-end vector and 



6 




200000 400000 600000 800000 1000000 

mcs 

FIG. 7; Simulation Results: End-to-end correlation function 
for component parallel to the walls for a chain with A'' — 
50, Rg = 53, and for slit widths of 20, 30, 40, and 5000 
(unconfined). 





Tube width (L) Simulation ti 


Predicted ti 


Ratio Sim/Pre 










N = 50 


20 


4.86E+06 


2.84E+05 


17.09 




30 


8.00E+05 


2.17E+05 


3.69 




40 


1.79E+05 


1.79E+05 


1.00 












N = 100 


50 


3.86E'i-06 


1.23E+06 


3.13 




70 


1.65E+06 


9.86E'i-05 


1.68 












N = 200 


80 


3.01 E+07 


7.22E+06 


4.17 




100 


1.78E+07 


6.22E+06 


2.86 



FIG. 8: Simulation Results: Observed value of longest odd 
mode versus predicted value. A factor of 16.75 is chosen as 
to match up the time scales at one point (A'^ = 50, L — 40). 
This factor absorbs any dependence on the monomer size b. 
The slit width is in units of lattice spacing. 
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(b) 

FIG. 9: Simulation Results: Time trace of the longitudinal 
component of the end-to-end vector P(t) for a chain with 
iV = 50, Rg = 53, and for slit widths of (a) 40 and (b) 30. 

time scales associated with individual monomer motion 
is predicted by the generalized Rouse dynamics, and seen 
in the simulations. 

In the next section, we analyze the motion of 
monomers using the framework of the Rouse modes. 

Monomer Motion 



the time scale increases as Rg. Under confinement, the 
slowest mode still involves reorientation, but now the 
time scale is expected to go as ~ N^L^^/"^, as ob- 
tained from our analysis of the Rouse modes. Based on 
these arguments, the reorientation events are expected to 
become increasingly rare as L is decreased. This expec- 
tation is borne out in our simulations, as evidenced by 
the time traces in Figl5] 

One way of describing fluctuations of the confined poly- 
mer is to say that the chain fluctuates about a fixed 
shape, characterized by the direction of P{t), chang- 
ing shape only occasionally. The fluctuations about 
the shape are described by the transverse Rouse modes, 
which as Eql2T] shows, relax faster as the slit width is 
decreased, and are independent of the size of the chain. 
Results from our simulations shown in Fig. [6] are con- 
sistent with these predictions. In simulations of a model 
polymer confined in a 2D box[2^, it has been observed 
that the slow modes correspond to shape changes, and 
that the slow dynamics is not reflected directly in the 
monomer motion. The shape of the polymer is a simpler 
variable in the slit geometry, but a distinct separation 
between time scales associated with shape changes and 



In this section, the effective Rouse modes are used 
to make detailed predictions for the mean-squared- 
displacemcnt (MSD) of monomers. The theory for the 
MSD of monomers should provide a framework for un- 
derstanding the dynamics inside biological cells where 
motion has been observed to be subdiffusive[27l| . 

Starting from Eq.® and making use of the fact the 
perpendicular mode index has a lower bound of N/ g, the 
number of blobs, it can be shown that the displacement 
squared of a monomer in the direction perpendicular to 
the walls is 



p=(N/g) 



1 — exp 



t 



(28) 



Likewise the expression for the parallel monomer motion 
is 



2t 4 

< (x„||(t)-x„||(0))2 >=j^+Y, 



p=i PI 



1 — exp 



t 



(29) 
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FIG. 10; Prediction for perpendicular motion from theory for 
a chain of 200 monomers, monomer size 2.8, and slit width 30 
{L/Rg = 1/5). (see Eq.(|28}) 




FIG. 11: Simulation Results: Plots of MSD of the perpen- 
dicular component of the middle monomer in a chain with 
= 200 The MSD have been rescaled by the maximum value 
of the saturation value at each L, and time has been scaled 
by, L"^"/"^, to capture the scaling of t^. These scalings lead 
to data collapse for L = 20 , 35 , 50. 

10OT rf 1 



The main features of the perpendicular motion, as pre- 
dicted by (^5)) , are depicted in Fig. pil)l . Since kp±_ has the 
same functional form as the unconfined polymer modes 
for p > N/g, the MSD exhibits a t^/^ power law that 
saturates to a value proportional to for t » r^, the 
relaxation time of the p = N/g mode. The explicit ex- 
ponent of 3/5 can be obtained by analyzing Eq. [28] for 
t <C T-'-. In this limit, the sums can be replaced by inte- 
grals, and we obtain: 




2l6/5i,4/3 

< {xn±{t) - x^AO)f >« 3^4/5 t"^T,{2/h) . (30) 



In the above equation, Tq is the Gamma function, and 

9 = 21/2^10/3 ■ The results of simulations for a iV = 200 
chain confined to slit widths ranging from 20 to 50 are 
shown in Fig. [TT] As predicted by the theory, MSD for 
different values of L can be collapsed on to each other by 
scaling the time t hy t-^ i^"/-^ and scaling the MSD 
by its maximum value which increases as (Fig. I12p . 

The MSD in the longitudinal direction, described by 
Eq. [211 is expected to exhibit multiple regimes. The 
functional form of of the effective spring constants, fcp||, 
changes at p = N/g with a increase at low p (for the 
odd modes) and a p^^"^ increase for the high p modes. 
The analysis is further complicated by the difference be- 
tween odd and even modes for low mode numbers. Even 
though there is no closed form expression for the MSD 
in the longitudinal direction, we can recover useful infor- 
mation about the functional forms and crossover times. 
To accomplish this, we look at times much less than the 
longest relaxation time of the system. With this assump- 
tion, the sum in (|29p can be approximated as an integral 



FIG. 12: Simulation Results: Plots of saturation value of 
displacement squared of the perpendicular component of the 
middle monomer in the chain recovered from simulations. The 
chain is 200 monomers long and plots are for slit widths of 
10, 20, 35, 50. 



and results in the following closed form expressions [2g 



< (a;„||(t) -x„||(0))2 > 

216/5^4/3 



2t 

N 



t'/'r ,4/3.., (2/5) + 



(31) 



27/2^2/3^1/3 



21/22,5/3 



where we have made explicit use of Eq. (|22|) . As ex- 
pected, the O*'' mode leads to purely diffusive motion. 
The system is dominated by the t^^^ behavior until a 
crossover time of associated with relaxations within the 
blob, tb\\: 



tb\\ 



2£m 

7r2fe4/3 



(32) 
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FIG. 13: Results for the MSD in the longitudinal direction 
obtained by numerical evaluation of Eq. (|29p for a chain of 200 
monomers, monomer size 2.8, and slit width 30. The dashed 
lines and the arrows mark the power laws and the crossover 
times predicted by the short-time approximation, Eq. (|32|) 



which is proportional to . For times longer than t^h, 
t^/"^ is the dominant behavior, and this subdifFusive mo- 
tion becomes diffusive as the Qth mode starts dominating 
at rc||[29,]: 



td\\ 



(33) 



These time scales, and the power laws are shown in 

Figdsi 

The different power-law regimes correspond to phys- 
ically different relaxation mechanisms. At the shortest 
time scales, the relaxation is not affected by the confine- 
ment, and the MSD follows that of the unconfined, self- 
avoiding chain. At intermediate times, the polymers can 
be thought of as an ideal chain of blobs, and blob com- 
pression dominates the dynamics [3]. The MSD, there- 
fore, shows the sub-diffusive behavior characteristic of 
an ideal chain with no excluded volume constraint. The 
long-time diffusive motion reflects the center-of-mass mo- 
tion. An interesting observation to be made is that the 
diffusive behavior sets in at earlier and earlier times as 
the slit width L is decreased. This apparently coun- 
terintuitive result, however, follows from the fact that 
with increasing confinement, transverse fluctuations de- 
cay faster, (r-*- ~ L^'^^^), and the longitudinal motion 
of the polymer essentially looks like that of a rigid rod. 
The reorientation events, which are the slowest dynam- 
ical process in the system, do not affect the MSD of 
monomers. 

The above analysis relies on our deductions for the 
functional form of the amplitudes which deter- 

mine fcp||. We have extracted the Rouse mode ampli- 
tudes from the simulations (see FigfT4|) to directly verify 
these functional forms. We find that for mode numbers 
p greater than a characteristic value Ps, the amplitudes 



decay as just as they do in the unconfined chain. 

However, for p < Ps, there is an observed splitting in the 
even and odd modes, with the even modes seemingly ex- 
hibiting a behavior and the odd modes decaying as 
Since we can obtain data only over a limited range 
(determined by the size of the chain), these power laws 
are only approximate. The crossover at Ps is a robust 
phenomenon. From our theoretical analysis, we expect 
Ps ^ N/g. For a chain of TV = 200, b 2.8 (character- 
istic of the bond- fluctuation model), and a slit width of 
L = 30, the blob model predicts that N/g = 8.0. The 
simulation data (Fig fM]) ) is in surprisingly good agree- 
ment with this prediction. Our simulation results, there- 
fore, support the conclusion that the statistical properties 
of the confined polymer can be understood in terms of 
the effective Rouse modes obtained from our theoretical 
framework. 
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FIG. 14: Simulation Results: Mode excitation vs mode num- 
ber p extracted from simulation for a chain 200 monomers 
long and a slit width of 30 lattice space. The saturation value 
for the largest modes is due to the chain being discrete instead 
of continuous. Solid and dashed lines represent the linearized 
Rouse theory predictions. 



We have used our simulations to carefully measure the 
MSD of monomers in the longitudinal direction, and the 
results are shown in Figs. [15] and 1161 Our results are 
consistent with the theoretical predictions. The short 
time behavior is governed by a sub-diffusive, t^^^ regime, 
associated with the unconfined polymer. There is a hint 
of an intermediate regime, governed by the slower power 
law of t^/^. Finally, at long times, the MSD is observed 
to increase linearly with t. The two characteristic times 
associated with the crossover 

fromi^/s to ^1/2 (tbii) and 
t^/^ to t (t£)||) are plotted in Fig [TBI and both times are 
observed to be monotonically decreasing with the tube 
widths, approximately obeying the scaling predicted by 
Eqs. (1221) and 
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FIG. 15: Simulation Results: MSD ol the middle monomer 
parallel to the slit walls for a chain of length 200 and the 
slit widths shown in the figure. For the unconfined chain, 
L = 5000, the diffusive regime is barely visible at t = 10^ 
mcs. In contrast, for L = 10, the diffusive regime is visible at 
t = lO'^ mcs. 



CONCLUSION 

We have presented a theoretical framework for ana- 
lyzing the dynamics of confined, self-avoiding polymers 
that is valid for any confining geometry. The formalism is 
based on renormalized, effective, Rouse modes. Through 
extensive simulations, we have shown that he predictions 
of the theory are borne out for a polymer chain confined 
to a slit. 

Given the Rouse mode framework, both collective dy- 
namics and single monomer motion can be determined 
from a knowledge of the mode relaxations. Our analysis 
shows that confinement introduces multiple time scales 
into the dynamics. The slowest mode is associated with 
the reorientation of the end-to-end vector, and its relax- 
ation time increases as L~^/^ with decreasing slit width. 
Transverse fluctuations of the chain relax on a time scale 
that decreases as 1/^°/^. This latter time scale shows up 
in the monomer dynamics as a crossover time from a t^^^ 
subdiffusive behavior, characteristic of the self-avoiding 
polymer, to a t^/^ subdiffusive behavior, characteristic of 
an ideal, phantom chain. In an unconfined polymer this 
intermediate regime is absent. Under confinement, the 
monomers exhibit pure diffusion in the longitudinal direc- 
tion at times longer than t^jh ~ L^^^, implying that the 
diffusive motion sets in at earlier times as the slit width is 
decreased. This observation can be understood by noting 
that a SAW behaves increasingly as a rigid rod as the slit 
width decreases and the dominant motion is that of diffu- 
sion of the center of mass. It should be remarked that in 
an unconfined chain the monomer motion becomes diffu- 
sive at time scales long compared to the reorientation, or 
rotational time scale. For the confined polymer, the sin- 
gle monomer motion does not exhibit the slowing down 
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FIG. 16: Simulation Results: (a) Plot of first crossover time 
from t^l^ to t^l^ (jB\\) from Figll5l(b) Plot of first crossover 
time from t^''^ to t (jd\\) from Fig |15l 



of the dynamics associated with reorientation events in 
which the chain undergoes a complete back bending, re- 
versing direction. The Rouse dynamics of an unconfined 
polymer is characterized by a single time scale, t\ ~ N'^" . 
As we have seen, the dynamics of a polymer confined to 
a slit has three distinct time scales, which scale with the 
slit width as L^°/^, L^/^, and L^"^!^. The separation be- 
tween the shortest, ^ L^^/"^ and the longest, L^'^l'^ time 
scales grows as L is decreased, resulting in a clear sepa- 
ration of the longitudinal and transverse dynamics. 

There are small discrepancies between theory and sim- 
ulations that can be ascribed to one or both of two ef- 
fects: (a) the lattice spacing associated with the sim- 
ulation model, (b) the finite length of the polymer, 
and the limitations of statistics. The one significant 
difference between theory and simulations is in the L- 
dependence of the longest relaxation time, t". Further 
studies are needed to investigate the origin of this discrep- 
ancy. Within the framework of Rouse modes, a more so- 
phisticated framework can be constructed by considering 
a renormalization of the single-monomer friction coeffi- 
cient, and its predictions compared to simulations. Ad- 
hering to the Rouse mode framework preserves the con- 
nection between static correlations and dynamics that is 
extremely desirable for parsing the complex dynamics of 
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polymers under confinement. 

The results presented in this paper demonstrate that 
the difference in dynamics in the longitudinal and trans- 
verse directions provides detailed information about the 
effect of tube confinement on a self-avoiding polymer. 
These results apply as long as the confinement is weak, 
i.e., the longitudinal dimension is larger than the uncon- 
fined radius of gyration of the polymer. If the longitu- 
dinal dimension is comparable or smaller than the ra- 
dius of gyration, then our results will be altered but the 
anisotropy, and especially the separation of time scales 
will persist and be prominent for aspect ratios larger than 
unity. The latter situation applies to bacterial chromo- 
somes. Recent research has focused on studying the role 
played by polymer confinement on chromosome dynam- 
ics and segregation foj. Our results indicate that resolving 
the transverse and longitudinal components in measure- 
ments of chromosomal dynamics would provide impor- 
tant diagnostic tools. Such measurements will aid the un- 
derstanding of the effects of cell geometry and crowding 
on the fine scale dynamics of bacterial chromosomes [2 7], 
and the relevance of intrinsic polymer dynamics to bio- 
logical processes. 

We would like to thank Jeremy Schmit and Jane Kon- 
dev for many insightful discussions, and Michael Hagan 
for a careful reading of the manuscript. This work was 
supported in part bv lnsf-dmr 04 03997 



[1] M. Daoud and J. P. Cotton, J. Phys. (Paris) 43, 531 
(1982). 

[2] F. Brochard and P.G. de Gennes, J. Chem. Phys. 67, 52 
(1977). 

[3] A. Aksimentiev, J. B. Heng, G. Timp and K. Schulten, 

Biophys. J. 87, 2086 (2004). 
[4] P.J. Lewis, Microbiology 147, 519 (2001). 
[5] K. Kremer, G.S. Grest, and I. Carmesin, Phys. Rev. Lett. 

61, 566 (1988). 



[6: 

[8 

[9 
[10 

[11 

[12 

[13 

[14 

[15 

[16 
[17 

[18 

[19 

[20 

[21 
[22; 

[23; 

[24; 

[25; 

[26; 

[27 

[28' 
[29 



G. Morrison and D. Thirumalai, J. Chem. Phys. 122, 
194907 (2005). 

R. M. Jendrejack, D. C. Schwartz, M. D. Graham and 

J.J. de Pablo, J. Chem. Phys. 119, 1165 (2003). 

K. Hagita and H. Takano, J. Phys. Soc. Japan 68, 401 

(1998). 

A. Arnold and S. Jun, Phys. Rev. E 76, 031901 (2007). 
M. Doi, Introduction to Polymer Physics (Oxford Uni- 
versity Press, 1995). 

A. Khoklov and A.I. Grosberg, Statistical Physics of 
Macromolecules (American Institute of Physics, 1993). 
K. Kremer and I. Carmesin, Macromolecules 21, 2819 
(1988). 

M. Doi and S.F. Edwards, The Theory of Polymer Dy- 
namics (Oxford University Press, 1988). 
D. Molin, A. Barbieri, and D. Leporini, J. Phys.iCond. 
Mat. 18, 7543 (2006). 

A. Arnold, B. Bozorgui, D. Frenkel, B. -Y. Ha and S. 

Jun, J. Chem. Phys. 127, 164904 (2007). 

P.J. Flory, J. Chem. Phys. 17, 303 (1948). 

S.F. Edwards and P. Singh, J. Chem. Soc. Faraday Trans. 

II 75, 1001 (1979). 

J.H. van Vliet and G. ten Brinke, J. Chem. Phys. 93, 
1436 (2005). 

R. Azuma and H. Takayama, cond-mat/9906021 (Jun 
1999). 

T. Sakaue and E. Raphael, Macromolecules 39, 2621 
(2006). 

P. Pincus, Macromolecules 9, 386 (1976). 

A. Di Cecca and J.J. Freire, Macromolecules 35, 2851 

(2002). 

C.-M. Chen and Y.-A. Fwu, Phys. Rev. E 63, 011506 

(2000) . 

N. Gilra, A.Z. Panagiotopoulos, and C. Cohen, J. Chem. 
Phys. 115, 1100 (2001). 

Y.-J. Sheng and M.-C. Wang, J. Chem. Phys. 114, 4724 

(2001) . 

A. Rahmanisisan, C. Castelnovo, J. Schmit, and C. Cha- 
mon (2007), aPS March Meeting. 

A. Fiebig, K. Keren and J. A. Theriot, Mol. Micro. 60, 
1164 (2006). 

J. Kalb, vol. Brandeis University (2008). 

P.G. de Gennes, J. Chem. Phys. 55, 572 (1971). 



